Consequences of the Pauli exclusion principle for the Bose-Einstein condensation of atoms and 

excitons 
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The bosonic atoms used in present day experiments on Bose-Einstein condensation are made up of fermionic 
electrons and nucleons. In this Letter we demonstrate how the Pauli exclusion principle for these constituents 
puts an upper limit on the Bose-Einstein-condensed fraction. Detailed numerical results are presented for hy- 
drogen atoms in a cubic volume and for excitons in semiconductors and semiconductor bilayer systems. The 
resulting condensate depletion scales differently from what one expects for bosons with a repulsive hard-core 
interaction. At high densities, Pauli exclusion results in significantly more condensate depletion. These results 
also shed a new light on the low condensed fraction in liquid helium II. 
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Recent experiments with ultracold fermionic gases have 
demonstrated the gradual crossover between a Bose-Einstein 
condensate of two-fermion molecules and a BCS-like conden- 
sate of fermion pairs 0113, Si- Turning this picture around, 
one might ask to what extent subatomic degrees of freedom 
play a role in Bose-Einstein condensates of bosonic atoms, 
because these atoms are made up of fermions: electrons and 
nucleons. From the energetic point of view there is no effect: 
subatomic excitation energies greatly exceed the thermal en- 
ergy scale of Bose-Einstein condensation. Therefore one can 
safely assume that the subatomic degrees of freedom are com- 
pletely frozen |4]. However, even for a frozen internal struc- 
ture one has to take into account the correct symmetries: at the 
level of the many-electron wave function, quantum mechan- 
ics dictates antisymmetry, which makes that electrons can not 
overlap. As a consequence, the Pauli principle for the elec- 
trons limits the available phase space for the bosonic atoms, 
which can have an influence on the properties of the conden- 
sate SHIlllllESIiil- It has been demonstrated before that 
a condensate of bosons made up of fermions has a maximum 
occupation number |12j]. For hydrogen atoms, that number 
corresponds to a condensate density of the order of 1/ (4tWq), 
with flo the Bohr radius. Such high densities are not reached in 
present day experiments on Bose-Einstein condensates 11311 . 
Still, the Pauli principle can have an effect also at lower den- 
sities, where it leads to condensate depletion. It is generally 
believed that it is sufficient to model this effect through an ef- 
fective interaction for the bosons which is strongly repulsive 
at short distances, like a hard-sphere potential or e.g. the (un- 
physical) r~ 12 term in the Lennard- Jones potential. The con- 
densate will be depleted, simply because of the excluded vol- 
ume. However, the only physical parameter which determines 
the low-density properties of the condensate is the scattering 
length. It is demonstrated below that any bosonic interaction 
with the right scattering length fails to reproduce the Pauli ex- 
clusion effect at high densities. We show how Pauli exclusion 
puts an upper bound on the Bose-Einstein condensed fraction 
of ultracold atomic gases. The bound is made quantitative for 
hydrogen atoms, through the use of an exactly solvable pair- 



ing model. The consequences for ultracold alkali gases, ex- 
citon condensates in semiconductors and liquid helium II are 
discussed. 

Following Penrose and Onsager fl7il . one can define a 
Bose-Einstein condensate by looking at the one-boson density 
matrix p£(r,r') of a many-boson system, 



p a (ry) = <* a |6j,6 r |*a), 



(1) 



with b T the operator that creates a boson at position r and <I>b 
the many-boson wave function. The system is said to exhibit 
Bose-Einstein condensation if the one-body density matrix 
has an eigenvector V|/fl(r) with an eigenvalue Xb of the same 
order as the total number of bosons, N. The ratio fg = Xb/N 
gives the condensed fraction, and the eigenvector VfB(r) cor- 
responds to the order parameter of the condensate. 

Taking into account that atomic bosons are actually made 
up of fermions, one realizes that the many-boson state |<J>g) 
corresponds at a more microscopic level to a many-fermion 
wave function For bosons made up of two fermions, the 
bosonic one-body density matrix of Eq.([0 can be related to 
the fermionic two-body density matrix. Grouping the fermion 
pair states in a single coordinate R = (ri,r2), one can write 
the fermionic two-body density matrix as a square matrix 
Pf (R,R'). A Bose-Einstein condensate would show up as an 
eigenvector \|/f(R) of the fermionic two-body density matrix 
with a macroscopic eigenvalue Xf |18]. If the bosons cor- 
respond to strongly bound pairs of fermions, one can expect 
the bosonic and the fermionic picture of the condensate to be 
equivalent |5], with 



\]f B (r)bldr = / \\t F (r u r 2 )al l al 1 dridr2, 



X R 



x f 



(2) 
(3) 



The fermionic model has the Pauli correlations between the 
fermions taken into account, while the bosonic model does 
not. Unfortunately, for any realistic model the fermionic 
many-body problem is too complicated in order to determine 
Pf (R,R') accurately. 
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A variational approach is feasible: X F is an eigenvalue of 
the many -body operator B^B, with the operator B^ defined as 



fit = / W(ri,r 2 )4 2 4 1 ^r 1 £/r 2 . 



(4) 



If one knows the structure of the order parameter of the 
fermionic pair condensate, one can determine an upper limit 
for the eigenvalue X F : 

X F = (<P F \B^B\<i> F ) <\E P \, (5) 

where Ep is the ground state energy of a fermionic pairing 
Hamiltonian, Hp = —B^B. This Hamiltonian is not meant to 
be phenomenological, but it is useful here because it is inte- 
grable fl9ll and exactly solvable j2(ill . Therefore, given the 
pair function \|/f(r), one can determine the ground state en- 
ergy Ep(2N) of Hp for 2N fermions, and make a rigorous 
variational statement about the boson condensed fraction: 

\E P {2N)\ 



N 



(6) 



This upper bound is valid at all temperatures, but will be most 
stringent at zero temperature, because there one expects the 
largest condensed fraction. The resulting eigenstate is not 
meant as an approximation to the true state \&f), because 
Eq.l|6j is a variational statement on the condensed fraction, 
not on the energy. 

To determine the upper bound of Eq.(|6}, one needs to know 
the bosonic order parameter in terms of the fermionic degrees 
of freedom. Although the following treatment is based on a 
hydrogen Is wave function, the procedure is more general. 
Alkali atoms such as 7 Li, 23 Na or 87 Rb can be seen as a paired 
state of a valence electron and a singly ionized core atom. Be- 
cause the tail of the wave function for the valence electron 
is similar to the hydrogen Is wave function, we expect qual- 
itatively the same Pauli effects as for hydrogen atoms. The 
Bloch-Messiah theorem says that it is always possible to write 
a fermion pair wave function in the form B^ = Li ¥f(O a [ a -i: 
where the index i does not necessarily refer to momentum 
states. The ground state energy of the corresponding pair- 
ing Hamiltonian Hp = B^B gives the upper bound for the 
condensed fraction. In the low-density limit this results in 
Sb < 1 -NZi\VF(i)\ 4 , while the upper bound in the high- 
density limit follows from the absolute maximum occupation 

number for the pair condensate 11211 : fg < 4^'^ ■ At in- 
termediate densities one can evaluate the exact solution for the 
ground state energy of Hp (see below). 

Here we consider hydrogen atoms in a cubic volume V with 
periodic boundaries, because for this case the wave function 
is known analytically, and therefore we can determine the up- 
per bound of Eq.Q easily. At low temperatures and densities, 
one can expect that the protons and electrons form hydrogen 
atoms and that all atoms are in a Is state. Because of transla- 
tional invariance, the bosonic condensate order parameter will 
be a uniform function in the center-of-mass coordinate of the 
atoms. The resulting pair operator can be written as 



(7) 



where the sum runs over all momentum states allowed by the 
periodic boundary conditions of the cubic volume V, with the 
subscripts p and e distinguishing between protons and elec- 
trons and \|/f(k) the pair wave function in momentum space, 
V|/f(k) = VZ/(1 +alk 2 ) 2 , with Z = 64%al/V. As discussed 
above, one can assume the internal structure of the hydrogen 
atoms to be completely frozen at temperatures of the order of 
1/jK, relevant for Bose-Einstein condensation |4]. Therefore 
we can rely on the pair operator of Eq.@ in order to evaluate 
the upper bound of Eq.(|6j. 

To find the ground-state energy of Hp for N pairs, one has 
to solve the following set of equations 1 2011: 
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with the function F(x) given by 
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(9) 



The energy Ep is given by Ep = Y,j=i -7 — 1 • The low density 
limit is obtained by taking the limit Z — > for a fixed number 
of pairs N. One obtains, for npa^ <C 1, that 



fs<l 



33% 



(n B al). 



(10) 



The result of Eq.( I10> was obtained without reference to in- 
teractions. At the bosonic level the Pauli blocking between 
the constituting fermions results in a repulsive interaction be- 
tween atoms at short distances. This interaction can be mod- 
eled using a pseudopotential 



V(r) 



4%h 2 a s d 
m or 



(11) 



with m the atomic mass. For non-interacting hard spheres, one 
can identify the scattering length a s with the radius lEH E2I1 . 
The Bogoliubov approximation, which applies in the low- 
density limit, results at zero temperature in a condensed frac- 
tion that scales as 122112311 



h = 1 - 



3 Vic 
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(12) 



The scaling with a power 1/2 assures that the variational 
bound Ea. (llO> is fulfilled. However, one observes that the Bo- 
goliubov result is fundamentally different from the expression 
of Ea. lllOt because it scales differently. This can be explained 
by the fact that the ground state wave function is not uniform. 
In fact, one can expect a higher amplitude for configurations 
where the bosons are well separated than for configurations 
where bosons nearly overlap 1 24] . Consequently, it turns out 
that Pauli blocking is only a second order effect in the low- 
density limit. The leading order is determined by bosonic 
many-body physics, which not only tries to avoid overlapping 
atoms, but also tries to minimize the energy. 
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The high density limit can be derived from the maximal 
occupation number of the pair state 11211 . which for hydro- 
gen results in fa < „ 1 » . At high densities one has to assess 

the effects of the interactions with the other atoms. Calcu- 
lations with a quantum Monte Carlo method l25ll and a con- 
strained variational method 1 26] have shown that for an inter- 
action with scattering length a s the condensed fraction is well 
described by the Bogoliubov approximation of Ea. (ll2t . up to 
densities of the order of ng ~ 10~ 3 a7 3 . For hydrogen atoms 
the s-wave scattering length is given by a s = 0.41ao 12711 . 
One can see in Fig. [2 that at a density of ng = 10~ 3 a7 3 = 
1.5 x 10~ 2 g!q 3 , the Pauli effect results in a condensate deple- 
tion of the order of 36%, while hard spheres with a radius 
a s , and hence any interaction with the same scattering length, 
yield a depletion of only 5%. One can conclude that at den- 
sities of the order of 10~ 3 «q 3 and higher the Pauli exclusion 
effect results in a significantly stronger condensate depletion 
than the hard-sphere potential, and that an effective two-body 
interaction for the bosons is not able to reproduce this effect. 

However, at high densities one also has to consider an- 
other effect: due to interactions with electrons and protons 
in neighboring atoms, the intrinsic wave function of the atoms 
might get deformed from the standard hydrogen Is wave func- 
tion. A way to take these interactions into account is to 
use a screened Coulomb potential, such as the Hulthen po- 
tential El: Vff(r) = - ^X-i ' with y o the hydrogen 
Is binding energy and y a dimensionless parameter propor- 
tional to the screening constant as defined in the Debeye- 
Hiickel or Thomas-Fermi models. This is a phenomenolog- 
ical way to take screening into account, which could de- 
viate from the true microscopic behavior at high densities. 
The intrinsic wave function in momentum space becomes 

W(k) = V^a 3 f)/[((l-Y) 2 +^)((l+Y) 2 +«^)]- 
In theory, at very high densities the system could undergo a 
Mott transition, where electrons and protons are no longer 
bound to each other in hydrogen atoms, but rather form a 
plasma. For cold hydrogen atoms this is not a realistic sce- 
nario, but it does apply to excitons in semiconductors, which 
have a similar intrinsic structure. The Mott transition occurs 
at an exciton density of n c ~ 0.02a,: 3 , with a x the exciton 
Bohr radius 1 30 ] . This behavior should be reflected in the 
density dependence of the screening parameter y: at very low 
densities y should tend to zero, in order to recover the hydro- 
gen wave function from the Hulthen wave function. The Mott 
transition, on the other hand, would require the bound states of 
the Hulthen potential to disappear around the transition den- 
sity. As this happens at y = 1, we follow Ref. I3H1 to take 
the screening parameter y = ^JriB/n c . The variational bound 
for the condensed fraction can be obtained for any density by 
solving the eigenvalue equations of the pairing Hamiltonian 
Hp with the screened pair structure. Fig.^shows the results 
for the hydrogen Is and Hulthen Is intrinsic wave functions. 
At the resolution of the figure the results for 1000 and 2000 
particles are indistinguishable, meaning that convergence has 
been reached and that finite-size effects are negligible at these 



particle numbers. The screening effects enhance the conden- 
sate depletion even more. They become important only when 
the maximal condensed fraction is lower than 80%, so they 
are of secondary importance in the high density regime. 
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FIG. 1: Maximal condensed fraction for hydrogen Is bosons (full 
line), Hulthen Is bosons (dashed line), and bosonic hard spheres of 
radius a s (circles), as a function of the density parameter iibOq- The 
hydrogen and Hulthen curves were calculated for 1000 bosons, the 
hard-sphere results were taken from Ref.|25]. 

The densities where the Pauli effect becomes sizeable are 
probably out of reach for an ultracold atomic hydrogen gas, 
because three-body recombination processes would convert 
the atoms into molecules. However, the same pair structure 
also applies to Wannier excitons in semiconductors. There 
the Pauli effect might explain, together with biexciton re- 
combination, why a clear signal of exciton condensation 
has not yet been observed in a three-dimensional structure. 
Bose-Einstein condensation of excitons has been observed 
in semiconducting bilave rs 1 1 411 . The physics there is ba- 
sically two-dimensional 1151 llal . There too the Pauli ef- 
fect applies. Although an analytical expression for the in- 
trinsic structure of these excitons is not readily available, 
one can estimate the Pauli effect by looking at the results 
of Fig. [2] for a two-dimensional hydrogen wave function, 
ff(t) °c (1 + a 2 k 2 /A)^l 2 and for a Gaussian wave function 
with the same low-density properties, « exp(— a 2 k 2 /5), 

with a x the two-dimensional excitonic Bohr radius. 

If one has to treat more valence electrons independently, 
then the boson operator corresponds to a three- or higher- 
body fermion operator instead of a pair operator, and the 
resulting Hamiltonian B^B is no longer exactly solvable. 
Still, one can expect the Pauli principle to have qualitatively 
the same effect on e.g. 4 He condensates: the electrons will 
avoid overlap and hence limit the available phase space for 
the 4 He atoms. A uniform distribution of hard-core bosons 
can qualitatively explain the reduced condensed fraction in 
liquid helium II compared to an ultracold low-density Bose 
gas 1171 12411 . The Pauli blocking of the underlying fermions 
offers a more microscopic view of this process. Pauli effects 
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FIG. 2: Maximal condensed fraction for two-dimensional excitons 
with a hydrogen-like (solid line) or Gaussian (dashed line) intrinsic 
wave function, as a function of the exciton density parameter n\,a\. 
The curves were calculated for 1000 excitons. 



will result in a significant depletion of the condensate at den- 
sities of the order of 10~ 3 times the close-packing density or 
higher, and definitely at the density of liquid helium. 

We have demonstrated here that Pauli blocking of the un- 
derlying electrons leads to condensate depletion in ultracold 
atomic gases and in exciton condensates. This effect might 
be measurable in systems where densities of the order of 
10 _3 aQ 3 can be reached, such as Wannier excitons in semi- 
conductors or 4 He films adsorbed on porous Vycor glass 13211 . 
This effect depends solely on the symmetry and internal struc- 
ture of the wave-functions. In the high-density regime, this 
effect can not be modeled through an effective two-body in- 
teraction at the bosonic level. Interactions might change the 
internal wave function of the fermionic pairs. Our calcula- 
tions based on a screened potential show that the Pauli effect 
dominates over the interaction effects at densities one or two 
orders of magnitude below the Mott transition density. 
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